Astronomy & Astrophysics manuscript no. 9058 


© ESO 2008 


February 5, 2008 





Chromospheric features of LQ Hydrae from Ha line profiles 

A. Frasca 1 , Zs. Kovari 2 , K. G. Strassmeier 3 , and K. Biazzo 1 

1 INAF - Catania Astrophysical Observatory, via S. Sofia 78, 1-95123, Catania, Italy 
e-mail: antonio.frasca@oact.inaf.it, katia.biazzo@oact.inaf.it 

2 Konkoly Observatory, H-1525 Budapest, P.O.Box 67, Hungary 
e-mail: kovari@konkoly . hu 

3 Astrophysical Institute Potsdam, An der Sternwarte 16, D- 14482 Potsdam, Germany 
e-mail: kstrassmeier@aip.de 

Received 14 Nov 2007 / Accepted 3 Jan 2008 

ABSTRACT 

We analyze the Ha spectral variability of the rapidly-rotating Kl -dwarf LQHya using high-resolution Ha spectra recorded during April-May 
2000. Chromospheric parameters were computed from the Ha profile as a function of rotational phase. We find that all these parameters vary 
in phase, with a higher chromospheric electron density coinciding with the maximum Ha emission. We find a clear rotational modulation of 
the Ha emission that is better emphasized by subtracting a reference photospheric template built up with a spectrum of a non-active star of 
the same spectral type. A geometrical plage model applied to the Ha variation curve allows us to derive the location of the active regions that 
come out to be close in longitude to the most pronounced photospheric spots found with Doppler imaging applied to the photospheric lines in 
the same spectra. Our analysis suggests that the Ha features observed in LQ Hya in 2000 are a scaled-up version of the solar plages as regards 
dimensions and/or flux contrast. No clear indication of chromospheric mass motions emerges. 
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1. Introduction 

According to the solar-stellar analogy, chromospheric stellar 
activity can be established by the presence of emission in the 
core of the Balmer Ha line. As in the Sun, Ha emission inten- 
sification has often been observed in surface features (plages) 
spatially connected with the photospheric star spots (see, e.g., 
Catalano et alil2000l: l5iazzo et al.ll2006l l2007l and references 
therein). Thus, the time variability of the Ha spectral features 
can be used to estimate the basic properties of the emitting 
sources, allowing the geometry of the stellar chromosphere to 
be mapped. 

In this paper we analyze the Ha spectral variability of 
the young, rapidly-rotating (P rot « 1.6 d) single K2-dwarf 
LQ Hydrae (HD 82558 = G1355) using 15 high-resolution 
spectra recorded during April-May 2000. All Ha spectra were 
acquired simultaneously with the mapping lines used for the 
year -2000 Doppler imag ing study presented in our previous pa- 
per dK6varietal.ll2004l hereafter Paper I). LQHya was first 
recognized as a chromospherically active s tar through Can 
H&K emission dBidelmanl 198 ll:lHeintzll 198 lh . The Ha absorp- 
tion filled i n by ch romospheric emission was reported first by 
Fekel et al] dl986h . and LQHya was classified as a BYDra- 
type spotted star. Variable Ha emission-peak asymmetry was 
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investigated by IStrassmeier et al. I JT993L hereafter Paper II), 
who attributed it to the presence of chromospheric velocity 
fields in the Ha forming layer probably surrounding photo- 
spheric spots. 

In our spectroscopic study in Paper I, we presented Doppler 
images using the Fe 1-6411, Fe 1-6430, and Ca 1-6439 lines 
for both the late April data and the early May 2000 dataset. 
Doppler imaging was supported by simultaneous photometric 
measurements in Johnson-Cousins VI bands. Doppler images 
showed spot activity uniformly at latitudes between -20° and 
+50°, sometimes with high-latitude appendages, but without 
a polar spot. Comparing the respective maps from two weeks 
apart, rapid spot evolution was detected, which was attributed 
to strong cross-talks between the neighboring surface features 
through magnetic reconnections. In Table[T]we give a summary 
of the stellar parameters as they emerged from Paper I. 

The observations are again briefly presented in Sect. [2] and 
the method used for the Ha spectral study is described in 
Sect. [3] The results are presented in Sect. H Since the Ha ob- 
servations of this paper were included in the spectroscopic data 
used in Paper I to reconstruct the year-2000 Doppler images, 
in Sect. [5] we take the opportunity to compare the photospheric 
features with the contemporaneous Ha-emitting regions. 
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Table 1. Astroph ysical data for LQHya (adopted from 
Kovari et alJl2004h 



Table 2. Observing log and radial velocities 



Parameter 


Value 


Classification 


K2 V 


Distance (Hipparcos) 


lo.jj+u.jj pc 


(B - V)ffipparcos 


0.933+ 0.021 mag 


(y — -Oflipparcos 


1.04+ 0.02 mag 


Luminosity, L 


0.270+0.009 L 




4.0+0.5 


Teff 


5070+100 K 


v sin / 


28.0+1.0 km/s 


Inclination, i 


65° + 10° 


Period, P rot 


1.60066+0.00013 days 


Radius, R 


0.97+0.07 R Q 


Mass 


« 0.8 M Q 


Age 


^ZAMS 



2. Spectroscopic observations 

The series of 15 Ha spectra was collected, one spectrum 
per night, during two observing runs (April 4-9, 2000, and 
April 25-May 3, 2000) at the Kitt Peak National Observatory 
(KPNO) with the 0.9 m coude-feed telescope. The 3096x1024 
F3KB CCD detector was employed, together with grating A, 
camera 5, and the long collimator. The spectra were centered 
at 6500 A with a wavelength range of 300 A. The effective res- 
olution was « 28 000 (11 km s" 1 ). We achieved signal-to-noise 
ratios (S/N) of about 250 in 45-min integration time, with the 
only exception being the spectrum at 2.726 phase for which an 
S/N of just 50 could be reached due to bad sky conditions. 

The mean HJDs of the observations, the phases, as well as 
the radial velocities, are summarized in Tabled For phasing the 
spectra, HJD = 2448270.0+ 1.600656 xE was used, where the 
arbitrarily chosen zero po int is the same as used for our former 
Dopp ler imaging studies dStrassmeier et al.|[l993l; iKovari et~aL 
l2004 . As radial velocity (RV) standard star, /3Gem (v r = 
3.23 km s" 1 ) was measured, except for HJD 245 1 643 .744 whe n 



16Vir (v r = 36.48 km s _1 ) was observed dScarfe et al.lll990h . 
For more details we refer to Paper I. 

3. Data analysis 

The B aimer Ha line is a useful and easily accessible indicator 
of chromospheric activity in the optical spectrum. It has been 
proven to be very effective for detecting chromospheric plages 
both in the Sun and in active stars, due to its high contrast 
against the surrounding chrom osphere (see, e.g.. iFrasca et al.l 
19981: iBiazzo et alJl2006Ll2007h . 

However, the Ha line is formed in a wide depth range in 
the stellar atmosphere, ranging from the temperature minimum 
to the upper chromosphere. To extract the chromospheric con- 
tribution from the line core, we apply a method usually called 
"spectral synthesis", which consists in the subtraction of syn- 
thetic spectra or of observed spe ctra of non-active standard 
stars (reference spectra) ( see. e.g..lHer big 1985; lBardenlll985l 



HJD 


Phase 


V r 


<T Vl 






(kms -1 ) 


(kms -1 ) 


2451639.681 


0.187 


7.1 


0.6 


2451640.647 


0.791 


7.8 


0.6 


2451641.648 


1.416 


8.6 


0.5 


2451642.662 


2.050 


7.4 


0.6 


2451643.744 


2.726 


7.6* 


2.4 


2451644.677 


3.309 


8.3 


0.6 


2451660.683 


13.308 


8.8 


0.5 


2451661.668 


13.924 


8.7 


0.6 


2451662.652 


14.538 


9.5 


0.6 


2451663.654 


15.165 


7.5 


0.5 


2451664.667 


15.798 


9.0 


0.6 


2451665.653 


16.413 


9.8 


0.5 


2451666.661 


17.043 


8.0 


0.5 


2451667.665 


17.670 


7.4 


0.6 


2451668.670 


18.298 


8.1 


0.7 



Frasca & Catalanolll994l: iMontes et al.l Il995h . The difference 
between observed and reference spectra provides the net Ha 



a with 16 Vir as RV standard 



emission, which can be integrated to estimate the total radia- 
tive chromospheric loss in the line. 

We point out that one can obtain the true chromospheric 
Ha emission with this method only if the chromosphere is also 
optically thin inside the active regions or if the plages are not 
extended enough to appreciably affect the photospheric spec- 
trum. In the opposite case, the subtraction of the underlying 
photospheric profile would overestimate the Ha chromospheric 
emission. 

For this reason we have also measured the Ha equivalent 
width in the observed spectra (EWn a ), fixing a constant in- 
tegration window of 5.0 A, i.e. integrating the emission core 
inside the absorption wings (see Fig. [T]). The net Ha equiva- 
lent width, AEWu a , was instead measured in the residual spec- 
tra obtained from the original ones with the "spectral synthe- 
sis" method. The non-active template was built using a high- 
resolution (R = 42 000) spectrum of the K2 V star HD 3765 re- 
trieved from the ELODIE archive. We chose H P 3765 beca use 
its spectral type and color index B - V = 0.942 dWeisll 1996b are 
nearly identical with those of LQ Hya , because its rotation ve - 
locity is very low (v sin i = 2.5 km s' MStrassmeier et aT]l2000|) . 
and because it has a low level of magnetic activity compared 
to LQ Hya. The spectrum was first degraded to the resolution 
of the KPNO spectra by convolving it with a Gaussian kernel 
and subsequently broadened by convolution with a rotational 
profile corresponding to the v sin i of 28 km s _1 of LQ Hya. 

The error of the equivalent width, cr EW , was evaluated by 
multiplying the integration range by the photometric error on 
each point. This was estimated by the standard deviation of the 
observed flux values on the difference spectra in two spectral 
regions near the Ha line. The values of EWn a , AEWn a , and 
(Tew are reported in Tabled 

We also measured the central wavelength of the Ha emis- 
sion both in the raw (observed) and residual spectra by evaluat- 
ing the centroid of the core emission. This quantity provided us 
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Fig. 1. Time series of Ha profiles for LQ Hya in April-May 
2000 sorted in phase order. The non-active template built up 
with a spectrum of the non-active, slowly-rotating K2V star 
HD 3765 is shown by a dash-dotted line at the bottom. The 
vertical dotted lines mark the integration interval for the raw 
Ha equivalent width (EWua)- 



with the radial velocity shifts between the Ha emission and the 
photosphere measured both in the raw (AV em ) and in the resid- 
ual spectra (AV^). The values of AV em and AV^ are reported 
in Table [3] as well. 



4. Results 

In the LQ Hya spectra, the Ha line is a strongly variable feature 
displaying absorption wings and a core changing from a com- 
pletely filled-in configuration to a moderate emission above 
the continuum (see Fig. Q]). A faint single peak just filling in 
the line core, without reaching the continuum, is observed at 
the phases of minimum Ha emission (0=0.2-0.4), whereas 
double-peaked emission profiles are observed with a central re- 
versal at the other phases. The "blue" emission peak is often 
stronger than the "red" one, although at some phases, nearly 
equal peak heights and even reversed intensity ratios are also 
observed. The Ha profiles are shown in Fig. [T] along with the 
photospheric template built up with an HD 3765 spectrum (see 
Sect. El). 



We estimated the chromospheric electron density in 
LQ Hya adopting the assumptions from Paper II, i.e. an isother- 
mal chromosphere with T c \x 0m = 10 000 K for which the opti- 
cal depth in the Ha cen ter can be derived by the formulation of 
Cram &Mullanl(l 1979k : 



/ A/ip ea k \ 2 

lnTchrom = i2A^j ' 



(1) 



where A/l peak is the wavelength separation of the blue and 
red emission peaks, and A/Id = 0.28 A is the chromospheric 
Doppler width. We measured AA pQ ^ on the observed spectra 
by fitting two Lorentzian functions to all the double-peaked 
Ha profiles. Obviously, this analysis could not be applied to the 
single-peaked spectra close to the minimum e mission at phase 
0.2-0 5. Combining equations 8 and 14 in ICram & Mullan 
the electron density n Q in cm 3 can be evaluated as 



n e = 1.67- 10 1 



^(^chrom) 



where is the ratio of the Ha peak to the continuum flux, 
as measured on our spectra, and B is the Planck function calcu- 
lated for the effective temperature of LQHya (r ef f = 5100K) 
and at r c hr m = 10 000 K. The optical depth at the Ha cen- 
ter changes from about 30 to 4 and the electron density from 
5-10 n to4-10 12 cm" 3 from phase ^0.2 (near the minimum Ha 
emission) to the phases of maximum emission (0=0.6-0.8, see 
Fig. [2]). In the same figure, the flux ratio between the red and the 
blue peaks |£ and the peak separation AA pe ^ are plotted ver- 
sus the rotational phase. All these quantities appear modulated 
by the stellar rotation. Symmetric profiles (|£ ^ 1) tend to be 
observed when the most active regions, judging from the pho- 
tospheric spot contrast, are in the visible hemisphere, similar to 
what was found by Strassmeier et al. in Paper II. 

The Ha equivalent width measured in the observed (EWu a ) 
and in the residual spectra (AEWna), the radial velocity shifts 
of the Ha emission (AVem and AV^), and the contempora- 
neous light curve in the V band (from Paper I) are plotted in 
Fig. [3] as functions of the rotational phase. A remarkable anti- 
correlation of V brightness and Ha emission is apparent. 

The highest electron density is observed when the Ha line 
is strongest, indicating a denser chromosphere above more ex- 
tended active regions, the opposite of what was observed in 
1991 (Paper II). We also found an average value for the chro- 
mospheric electron density that was higher than the one re- 
ported in Paper II (2 • 10 11 cm" 3 ). Moreover, the rotational mod- 
ulation of all the quantities deduced by the Ha line analysis is 
much more evident in the present data. Presumably, in April- 
May 2000, LQHya was more active than in 1991, and it dis- 
played more and larger active regions. 

5. Discussion and conclusions 

To obtain information about the surface location of the ac- 
tive regions in the chromosphere of LQ Hya, we applied a 
simple geometric plage model to the rotationally modulated 
chromospheric emission. This method had been described and 
successfully applied to the Ha modulation curves of sev- 
eral other active stars (e.g., iFrasca et al. I l2000l iBiazzo et al 
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Fig. 2. From top to bottom. The computed electron density n Q 
(cm -3 ), the flux ratio between the red and the blue peaks, and 
the peak separation in A against the rotational phase. 



20061: iFrasca et all 12007k . On LQ Hya, two circular bright 



spots (plages) are normally sufficient for r eproducing the ob 



served variations withi n the data errors (cf . F rasca et al. hood 



iKovari & Bartusll2001[) . In our model the flux ratio between 
plages and the surrounding chromosphere (^piage/^chrom) is 
a free parameter. Solar values of ^piage/^chrom ~ 2 , as de- 



duced from averaging many plages in Ha (e.g., lEUisonl 



19521: lLaBonte1ll986l:lAvres et al.lll986h . are too low to model 
the high amplitudes of Ha AEW cur ves in very active stars 
dBiazzo et al.ll2006l : lFrasca et~aT1l2007 ). In fact, extremely large 
plages, covering a significant fraction of the stellar surface, 
would be required with such a low flux ratio, and they could 
not reproduce the observed modulation. In order to achieve a 
good fit, a flux ratio of 5, which is also typical of the brightest 



Svestk; 



>ngntest 

2tS; 



parts of so lar plages or of flare regions (e.g., 
IZirinlll988h . was adopted. 

The solutions essentially provide the longitude of the 
plages and give only rough estimates of their latitude and size. 
The information about the latitude of surface features can be 
recovered from the analysis of spectral-line profiles broadened 
by the stellar rotation, as we did for the spots of LQ Hya with 
Doppler imaging in Paper I. A similar technique cannot reach 
the same level of accuracy when applied to the Ha profile 
whose broadening is dominated by chromospheric heating ef- 
fects, which are particu larly efficient in the active regions (e.g., 
Lanzafame et al.ll2000h . 



We searched for the best solution by varying the longi- 
tudes, latitudes, and radii of the active regions. The radii are 



Fig. 3. Top panel: V light curve contemporaneous to the spec- 
tra. Middle panels: Rotational modulation of the raw Ha equiv- 
alent width (EWua, open circles) and of the residual Ha EW 
(AEWu a , filled circles). Bottom panel: Velocity shift between 
Ha emission and photospheric lines both for observed (open 
circles) and residual spectra (filled circles). The continuous 
lines in the middle and bottom panels represent the best fit of 
our 3 -plage model (cf. Fig. [4]), while the dotted lines represent 
the result of our 2-plage model. 

strongly dependent on the assumed flux contrast ^piage/^chrom- 
Thus, only the combined information between plage dimen- 
sions and flux contrast, i.e. some kind of plage "luminosity" 
in units of the quiet chromosphere (Lpiage/Ajuiet)* can be de- 
duced as a meaningful parameter. Note also that we cannot 
estimate the true quiet chromospheric contribution (network), 
since the Ha minimum value, A2sWq U iet = 0.84 A, could still 
be affected by a homogeneous distribution of smaller plages. 
However, such an approximation seems valid for LQ Hya be- 
cause our aim was only to compare the spatial distribution of 
the main surface inhomogeneities at chromospheric and photo- 
spheric levels. 

We also searched for solutions with three active regions, 
finding a small but possibly significant increase in the good- 
ness of the fit, with the x 2 passing from 5.77 to 4.45 (see 
Fig. [3]). The model with only two plages requires features at 
rather high latitudes (Table 0]) to explain the nearly flat and 
long-lasting maximum, while a 3 -plage model allows placing 
plages at lower latitudes, in better agreement with the spot lo- 
cations from Doppler imaging. However, the longitudes of the 
two largest plages do not change very much (less than 15 de- 
grees). 
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Table 3. Equivalent widths of the raw (EWn a ) and the resid- 
ual spectra (AEWu a ), errors of the equivalent widths (<tewX 
and velocity shifts AV em and A for the raw and the residual 
spectra, respectively. 



Table 4. Plage parameters for LQ Hya in April-May 2000. 



Phase EW Ua AEW Ua 0~ EW 





(A) 


(A) 


(A) 


(kms- 1 ) 


(kms" 1 ) 


0.187 


0.172 


0.934 


0.041 


-27.46 


-1.78 


0.791 


0.117 


1.003 


0.038 


-14.05 


-0.13 


1.416 


0.183 


0.935 


0.036 


-29.11 


-0.96 


2.050 


0.107 


1.013 


0.043 


7.40 


-1.29 


2.726 


0.116 


1.045 


0.213 


-17.00 


-5.80 


3.309 


0.239 


0.879 


0.037 


-32.70 


-2.92 


13.308 


0.248 


0.858 


0.044 


-28.68 


-2.19 


13.924 


0.058 


1.045 


0.048 


-10.77 


1.38 


14.538 


0.025 


1.050 


0.041 


-20.01 


-1.38 


15.165 


0.099 


0.998 


0.041 


-0.69 


-0.56 


15.798 


0.073 


1.050 


0.036 


-1.87 


1.09 


16.413 


0.215 


0.919 


0.035 


-29.60 


-0.49 


17.043 


0.094 


1.024 


0.034 


-0.18 


0.65 


17.670 


0.038 


1.074 


0.051 


-13.38 


0.36 


18.298 


0.325 


0.834 


0.038 


-43.07 


0.29 



The resulting plage longitudes are in good agreement with 
the Doppler results in Paper I; i.e., photospheric minima at 0.6 
and 0.9 with the largest/coolest photospheric spots correspond 
to the most luminous chromospheric phases, and vice versa, 
the brightest photometric phase at 0.3 overlap with the low- 
est chromospheric Ha emission, again supporting the paradigm 
that photospheric spots are physically connected with chromo- 
spheric plages (Fig. [4]). 

The synthetic Ha EW curve for the model with three plages 
is plotted in the middle panel of Fig. [3] as a continuous line su- 
perimposed on the data (dots), while the 2-plage solution is 
displayed by a dotted line. Both of them reproduce the ob- 
served rotational modulation quite well. With the same models 
we were able to calculate the radial velocity shift between the 
Ha synthetic emission profile, resulting from the quiet chromo- 
sphere plus the visible plages, and the photosphere. It is clear in 
the bottom panel of Fig. [3] that the amplitude of the theoretical 
AVem curves, both for a 2-plage and a 3-plage model, is very 
low, consistent with the values derived from the residual pro- 
files (AVJ^) and in total disagreement with the velocity shifts 
derived from the raw spectra. This strongly supports the spec- 
tral synthesis method we used to evaluate the chromospheric 
emission in the Ha line. 

We would like to outline the rotational modulation of other 
features observed in the Ha profile, such as the peak separa- 
tion, which is related to the chromospheric electron density 
(Sect.H]) and the asymmetry of blue/red peak intensity (Fig. [2]). 
The chromosphere of LQ Hya displays a higher electron den- 
sity above active regions. The blue emission peak is stronger 
than the red one (|^ < 1) in the "quiet" chromosphere of 
LQ Hya, while nearly equal peaks (|£ ^ 1) tend to be observed 
when the most active regions are in the visible hemisphere. 



Radius Lon. Lat. 



^chrom 



(°) 


O (°) 




2 Plages Solution 


15.5 


215 35 


5 


20.0 


353 66 


5 


3 Plages Solution 


13.0 


210 20 


5 


13.0 


8 15 


5 


9.5 


280 15 


5 



Asymmetry in the peaks of an emission line with a central 
reversal has b een frequently observed for Can K l ine both in 
the Sun (e.g..lzirin|[l988l: iDing & Schleicherll 19981 and refer- 
ence therein) and in cool stars (e.g.. IMontes et al 2000, and 
reference therein). The blue asymmetry is more frequently ob- 
served in the solar chromosphere and is commonly a ttributed 
to the propagation of acoustic waves (e.g.. lCramlll976l) with an 
upward velocity on the order of 10 km s" 1 in the layer in which 
the K2 emission peaks are formed or a downw ard motion in 
the la yer producing the central reversal K3 (e.g. JPurrant et al 
1197 ' 



ll926|). 
lOranie 



lOranjel (11983b has shown that the average Ca 11 K emission 
profile in solar plages has a completely different shape from the 
surrounding chromosphere, with the red peak slightly stronger 
than the blue one. This could reflect the different physical con- 
ditions and velocity fields in active regions compared to the 
quiet chromosphere. A similar analysis cannot be made for the 
Sun as a star in Ha, because this line is an absorption feature 
in the quiet chromosphere and only a filling in is observed in 
plages. However, similar changes of the Ha line profile can be 
expected in the plages of any stars that are more active than the 
Sun. Thus, the rotational modulation of ^ is consistent with 
the presence of plages in the chromosphere of LQ Hya. 

These results suggest scaled-up versions (regarding size 
and brightness) of solar-type plages for the Ha features ob- 
served in LQ Hya in 2000, without any clear evidence of strong 
mass motions in its chromosphere. 
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Fig. 4. Schematic representation of the surface features of 
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